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We provide a "Baedecker" or travel guide to the directions on the sky where the dark matter 
annihilation signal may be expected. We calculate the flux of high energy 7-rays from annihilation 
of neutralino dark matter in the centre of the Milky Way and the three nearest dwarf spheroidals 
(Sagittarius, Draco and Canis Major), using realistic models of the dark matter distribution. Other 
investigators have used cusped dark halo profiles (such as the Navarro-Frenk- White) to claim a 
significant signal. This ignores the substantial astrophysical evidence that the Milky Way is not 
dark-matter dominated in the inner regions. We show that the annihilation signal from the Galactic 
Centre falls by two orders of magnitude on substituting a cored dark matter density profile for 
■ a cusped one. The present and future generation of high energy 7-ray detectors, whether atmo- 

' spheric Cerenkov telescopes or space missions like GLAST, lack the sensitivity to detect any of the 

monochromatic 7-ray annihilation lines. The continuum 7-ray signal above 1 GeV and above 50 
GeV may however be detectable either from the dwarf spheroidals or from the Milky Way itself. If 
5— ( ' the density profiles of the dwarf spheroidals are cusped, then the best prospects are for detecting 

Sagittarius and Canis Major. However, if the dwarf spheroidals have milder, cored profiles, then 
the annihilation signal is not detectable. For GLAST, an attractive strategy is to exploit the wide 
field of view and observe the Milky Way at medium latitudes, as suggested by Stoehr et al. This is 
reasonably robust against changes in the density profile. 



I. INTRODUCTION 
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■"sj" . The foremost candidate for the cold dark matter (CDM) composing galactic haloes is the lightest neutral supersym- 
metric particle, namely the neutralino 1]. If so, then neutralino pair annihilation may lead to observable consequences, 
in particular the emission of high energy 7-radiation 2J . The possibility that such 7-rays may be identified by forth- 
coming atmospheric Cerenkov telescopesfACT) such as VERITAS (3| or by satellite-borne detectors like GLAST Q 
] has excited considerable recent interest IE S 13 • 

■ It is clearly of importance to identify the best places to search for such an annihilation signal. Inspired by the 
highly cusped models based on numerical simulations of dark halo formation |^, Il0| , a number of investigators have 
II suggested that the centre of the Milky Way may be the optimum target. For example, Bergstrom et al. have 
Q shown that if the dark matter density is cusped as l/r at small radii, then the 7-ray flux would be detectable for 
^ , typical neutralino properties in the minimal supersymmetric extension of the Standard Model. Inspired also by the 
^ • persistence of substructure in numerical simulations, a number of authors ^5,, 11, 12j have argued that a substantial 
\ enhancement in the 7-ray signal can be expected from such 'clumps'. In these calculations, the inner regions of the 
^ ' substructure are also usually assumed to be cusped. However, even within the framework of the cusped models favored 
by cosmological simulations, these conclusions have been contested as being overly optimistic 6]. 

More awkwardly, there is a substantial body of astrophysical evidence that the halo of the Milky Way is not cusped 
c3 ■ at all [3. First, the microlensing optical depth towards the Galactic Genter is very high. Particle dark matter does 
not cause microlensing, whereas faint stars and brown dwarfs do. The total amount of all matter within the Solar circle 
is constrained by the rotation curve, so this tells us that lines of sight towards the Galactic Genter are not dominated 
by particle dark matter. More specifically, haloes as strongly cusped as l/r, normalised to the local dark matter 
density as inferred from the stellar kinematics in the solar neighbourhood, are ruled out by the high microlensing 
optical depth 14]. Second, the pattern speed of the Galactic bar is known to be fast from hydrodynamical modelling 
of the motions of neutral and ionised gas. If dark matter dominates the central regions of the Milky Way, then 
dynamical friction will strongly couple the dark matter to the Galactic bar and cause it to decelerate on a few bar 
rotation timescales TSI. It is now largely accepted by astronomers that bright galaxies like the Milky Way do not 
have cusped dark haloes today, with some investigators suggesting that feedback from star formation may provide a 
resolution with cold dark matter theories 0| . 

In fact, there is no observational evidence whatsoever that any nearby galaxy has a cusped dark halo profile. The 
rotation curves of low surface brightness and dwarf spiral galaxies have been the subject of a long controversy [TtLITsI I. 
The effects of beam smearing mean that the H I rotation curves of many dwarf spirals are broadly compatible with 
both cores and cusps. However, the H II rotation curves for at least some dwarf spirals are not compatible with cusps 
[T9I I . Most dwarf spheroidals (dSphs) do not contain gas and so the structure of the dark halos must be inferred from 
stellar motions. Very recently, the survival of kinematically cold substructure in the Ursa Minor dSph has been used 
to argue against a cusped halo (221]. Hence, even at the least massive and most dark matter dominated end of the 
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galaxy mass spectrum, the predictions of cold dark matter theories concerning halo structure seem to disagree with 
the observations. Nonetheless, for the three nearest dSphs - Draco, Sagittarius and Canis Major - there is no direct 
evidence either for or against central cusps in their dark matter distribution. 

Given this weight of evidence, it seems very prudent to use both cored and cusped halo models to estimate the 
range of the expected 7-ray annihilation signal. In this paper, we examine four possible locations - the Galaxy Center 
and the centers of the three nearest dark-matter dominated dwarf spheroidals (Draco, Sagittarius and Canis Major). 
In Section II we use the most recent data on the velocity dispersion of the dSphs to constrain a variety of dark halo 
models, and in Section III evaluate the 7-ray flux from neutralino self-annihilations. Sections IV and V summarise 
the expected contribution from the background and the criterion for detection respectively. The results for second 
generation ACTs and for GLAST are given in Section VI. 



II. MODELS OF DWARF SPHEROIDALS 



Dwarf spheroidals (dSphs) warrant attention because they are amongst the most extreme dark matter dominated 
environments. For example, the mass-to-light ratio of Draco is ~ 250 in Solar units [23|, while that of the Sagittarius 
is ~ 100 22]. The recently discovered possible dSph in Canis Major seems similar to the Sagittarius in structural 
properties and dark matter content '23l. Given the seeming absence of dark matter in globular clusters, dSphs are 
also the smallest systems dominated by dark matter. 

We develop two sets of models of dSphs. The first set is the cored spherical power-law models p^ : 

2 Q 3 2 I 2Q _ \ 
Ppowlr-j- (^2+^2)2+a/2 ' 

Here, Vc is the core radius and Va is a velocity scale. When a = 0, the model has an asymptotically flat rotation curve 
and is the cored isothermal sphere. The rotation curves of dwarf galaxies may be gently rising or falling at large radii, 
so we also consider models with a — —0.2 and 0.2 respectively. 
The second set of models is the cusped haloes 

Pcusp = — ^- rrzz, (2) 

favored by numerical simulations. Here, is the scale radius and A is the overall normalisation. When 7 = 1.5, the 
model is the highly cusped Moore et al. profile, when 7 = 1, the model is the Navarro- Frenk- White (NFW) profile 
[To| . Additionally, we study the case 7 = 0.5 which represents a still milder cusped profile. 

The two free parameters determining the shape of the profile are set by fitting to observational data on the Draco 
dSph using the Jeans equation l25l|. For a spherical galaxy, the enclosed mass M{r) is related to observables via 

, , rivr)"^ ( dlogi^ dlog(t;,^) \ , , 

G \dlogr dlogr / 

Here, v is the luminosity density, (w^) is the radial velocity dispersion of the stars and /3 is the anisotropy of the stellar 
motions. The luminosity density of Draco v is taken as |2a |: 

-- (4) 



(rg+r2)5/2' 



with ro = 9.71' « 0.23 kpc (using a heliocentric distance for Draco of 82 kpc). There are 6 observational points showing 
the line of sight velocity dispersion of Draco at different radii in |23| (using the data with no rotation subtracted). The 
datapoints are consistent with a flat profile between 2' and 22'. Assuming that the anisotropy now vanishes, then the 
radial velocity dispersion is equal to the linc-of-sight velocity dispersion. Finally, the left-hand side of eq. is fitted 
to the known right-hand side at the locations of the datapoints in J7], thus giving estimates for the two unknown 
parameters for each density profile as quoted in Tabled 

This algorithm provides models of the Draco dSph that satisfy the available observational data. Unfortunately, the 
radial variation of the velocity dispersion has not been measured for the Sagittarius dSph. However, the central line- 
of-sight velocity dispersion of Sagittarius dSph is 11.4 kms"""^, very similar to that of Draco (10 kms""'^). Henceforth, 
we assume that the underlying structural parameters {vz,-,Tc or A, rg) of the Sagittarius dSph are the same as Draco. 
The third dSph under study - Canis Major - has only recently been claimed and the evidence for its existence is not 
yet clear-cut. There is certainly a surprising concentration of stars in the direction of Canis Major, but this could 
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Cored Power-Law Models 



a 


Va km s ^ 


Tc kpc 


rt kpc 
MW - Iso 


rt kpc 
MW - NFW 


M(rt) ~ 1O«M0 


0.2 


24.7 


0.25 


6.2 (2.16) 


1.3 (0.5) 


4.6 (2.0) 





22.9 


0.23 


7.8 (2.5) 


1.4 (0.51) 


9.5 (3.0) 


-0.2 


20.9 


0.21 


10.1 (2.8) 


1.6 (0.52) 


22.43 (4.9) 



Cusped Models 



7 


A X IO'^Mq 


rs kpc 


n kpc 


rt kpc 


M{rt) ^ W^Mq 








MW - Iso 


MW - NFW 




0.5 


2.3 


0.32 


6.6 (2.5) 


1.5 (0.6) 


5.5 (3.1) 


1 (NFW) 


3.3 


0.62 


7.0 (2.59) 


1.6 (0.57) 


6.6 (3.5) 


1.5 (Moore) 


2.9 


1. 


6.5 (2.4) 


1.5 (0.6) 


5.5 (2.8) 



TABLE I: Parameters of the dark matter halo profiles of the Draco dSph. The last three columns also give (in parentheses) 
the values at the location of the Sagittarius dSph. Two values are given for the tidal radius, according to whether the Milky 
Way halo is modelled with an isothermal power-law model or a NFW model. [Notes: (1) the models sometimes require a slight 
velocity anisotropy in the very innermost parts to ensure everywhere physical stresses in the Jeans equation, (2) the scale radius 
rs is constrained to lie below 1 kpc]. 



be due to some dynamical feature like an outer spiral arm associated with the Milky Way disk. However, for the 
purposes of this study, we assume the interpretation of the data in terms of a merging dSph is correct and that Canis 
Major dSph is similar in size and structure to Sagittarius. 

To determine the extent of the dark matter halo of the dSphs, the tidal radius must be estimated. The approximate 
method used conventionally is derived from the Roche criterion. The tidal radius is found by requiring that the 
average mass in the dSph is equal to the average interior mass in the Milky Way halo, namely 

^dsph(?'t) _ MMw(?'dSph - n) 
rt (?'dSph ~rt)^ 

Here, Muwir) and M<iSph('') are the masses enclosed within radius r of the Milky Way halo and the dwarf spheroidal 
respectively, while rdSph is the distance from the Galactic Center to the centre of the dSph. We remark that this is 
not the same as the procedure used in a number of recent papers 0, , in which the local density at the tidal 

radius in the dSph is set equal to the density of the background Milky Way halo at the center of the dSph. 

The results depend on the choice of profile for the Milky Way halo. For comparison purposes, we consider both 
a cored isothermal profile with rc = 10 kpc and Va = 220 kms~^, and a NFW profile with concentration parameter 
c = 10. The total mass of the Milky Way halo is fixed at Mmw ^ 10^^ A^©, as suggested in |33|. We show in Tabled 
the results when eq. Q is used to determine the tidal radius of a dSph at the locations of Draco and Sagittarius, for 
the two adopted models of the Milky Way halo. 



III. THE GAMMA-RAY FLUX 



Let the neutralino mass be and its self-annihilation cross-section be (crw). Then, the 7-ray flux from neutralino 
annihilation is given by 

= X -1- / dnf pMs)] ds, (6) 

47rm^ An J^n Jios 

where p is the density of the dSph as a function of distance from its center r, which of course depends on the 
heliocentric distance s. The integration is performed along the line-of-sight to the target and averaged over the solid 
angle Ail of the detector. In particular, — 2 for the annihilation of two non-relativistic neutralinos into two 
photons (xx ^ 77) and = 1 for the annihilation into a photon and a Z boson (xx ~* ^l)- The first part of 
the integrand © depends on the particular particle physics model for neutralino annihilations. The second part is 
a line-of-sight integration through the dark matter density distribution. We discuss each in detail in the next two 
subsections. 
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mSUGRA parameters 



mo (GeV) 


mi/2 (GeV) 


tan P 


1^0 1 (GeV) 


10-10000 


10-10000 


1-60 


10-10000 



TABLE II: The portion of the mSUGRA parameter space randomly scanned to generate the models. Here, mo and mi/2 are 
respectively the common scalar and gaugino mass at the unification scale, while Ao is the trilinear parameter and tan /3 is the 
ratio of the vacuum expectation values of the two Higgs fields. The term in the Lagrangian is allowed to have either sign. 



A. Particle Physics Model 

To compute N^{av) j [A'Km'^, we have to select a supersymmetric model. We focus on minimal supergravity 
(mSUGRA) models with universal gaug ino and scalar masses and trilinear terms at the unification scale jS l| . We use 
the computer programme SoftSusy |32l| to scan the supersymmetric parameter space (see Table and generate 10^ 
models which have consistent electroweak symmetry breaking and grand unification. The output at the electroweak 
scale is fed into the programme DarkSusy which computes the relic density and products of the neutralino 
annihilations. It also checks that a given model is not ruled out by present accelerator experiments. 

A feasible model is one which is permitted by accelerator limits and which predicts a relic density in the range 
0.005 < ficDM^^ < 0.2. This is somewhat broader than the range 0.09 < ricDM^i^ < 0.13 determined by fitting the 
standard ACDM model to the WMAP data ,341 . This is done so as to incorporate the higher values for f^cDM^^ 
found for consistent alternative CDM models j^. The lower limit is set by requiring the relic particle to provide 
most of the dark matter in galaxies (taking a typical mass-to-light ratio for galaxies of ^ 10, cf. the critical mass to 
light ratio of ^ 2000 in solar units). The range of the supersymmetric parameters scanned are given in Table ITU For 
each feasible model, we record the quantities N-y{av) for the discrete lines XX ~^-J^ ^^'^ XX ~^ ■^7) ^.s well as the 
continuous 7-ray spectrum above 1 and 50 GeV. Note that in the previous study , the relic density was much less 
constrained and models with an arbitrary low values were permitted. 



B. Line-of-Sight Integration of Dark Matter Density 

The line-of-sight integration can be manipulated thus: 



(J) AO = Xo / JW'^^ = 1^ / '^^ '^""^ / """^ '^'^P^ Us^ + sl~2sso cosd) (7) 

where 

JW^ I dsp\r). (8) 

J los 

In these formulae, angled brackets denote the averaging over the solid angle Afi, while Smin ^^nd Smax are the lower and 



upper limits of the line-of-sight integration, given by sq cos 9 ± y — Sq sin 6. Here, sq is the heliocentric distance 
of the dSph and rt is the tidal radius of the dSph. Finally, 0„iax is the angle over which we average around the center 
of the dSph. It generally is, at least, equal to the experimental resolution and can be fixed using: 



An^2n [ de sine = 27r(l - cos(e'max))- 
^0 



The quoted point spread function widths for the various experiments are: 0.4° (EGRET), 0.1° (GLAST, HESS and 
VERITAS), 0.15° - 0.04° (CANGAROO-III). EGRET and GLAST are sateUite detectors with low energy thresholds 
(« 100 MeV), high energy resolution (w 15%) but only moderate angular precision. The others are ACTs with higher 
thresholds (« 100 GeV) but better angular resolution. Typical reference sizes for the solid angle are Ail = lO"'^ sr 
for ACTs and GLAST and An = IQ-'^ sr for EGRET. 

Table HTll hsts values of (J) An for the dSph profiles introduced in Section XL The heliocentric distances to the Draco, 
Sagittarius and Canis Major dSphs are ~ 80, 24 and 8 kpc respectively and this largely controls the relative values 
of (J)Aa for the three dSphs. Clearly, the comparative closeness of the Canis Major dSph works to its advantage as 
a possible target. An increase in angular sensitivity enhances the signal for all three dSphs. We remark that, in the 
literature, there is a considerable spread in the values obtained for (J)Aa for different sources. In 0, a 1-component 
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Cored Power-Law Models 



a 


Sagittarius 


Draco 


Canis 


Afi = 10"^ sr 


An = 10"^ sr 


Af7 = 10"^ sr 


Af2 = 10"^ sr 


Afi = 10"^ sr 


Af7 = 10"^ sr 


0.2 


0.6 


3.4 


0.07 


2.2 


2.4 


3.4 





0.6 


3.3 


0.06 


2.2 


2.4 


3.5 


-0.2 


0.6 


3.2 


0.07 


2.2 


2.4 


3.4 



Cusped Models 



7 


Sagittarius 


Draco 


Canis 


An = 10"^ sr 


AVL = 10"^ sr 


Afi = 10"^ sr 


Afi = 10"^ sr 


AO = 10"^ sr 


AO = 10"^ sr 


0.5 


1.1 


17.8 


0.1 


5.7 


6.2 


32.3 


1 (NFW) 


1.3 


36.9 


0.1 


7.2 


8.3 


139.9 


1.5 (Moore) 


7.3 


615.1 


0.6 


55.4 


49.1 


5469 



Galactic Center 



Profile 


AO = 10"^ sr 


AO = IQ-'^ sr 


NFW, 7 = 1 


26 


280 


Cored, a = 


0.3 


0.3 



TABLE IIL Values of (J) An for the Sagittarius, Draco and Canis Major dSphs in units of 10^^ GeV^cm"^. The tidal radius 
of the dSphs is calculated assuming an isothermal profile for the Galactic halo. Additionally, results in the direction of the 
Galactic Center are given for both the NFW and isothermal models of the Galactic halo. 



King profile was used to model the dSph density distribution. These authors only give explicit estimates of the entire 
linc-of-sight integral. However, the values are of the order of lO^^GeV^cm"^, lower than those implied by Table ITTTI 
In [8|, no angular average is taken, but instead the approximation 

is used, together with an singular isothermal profile for Draco. Taking the value given for — 100 GeV implies 
3.7 X 10^" GeV^ cm~^ for the line-of-sight integral. It is surprisingly low and yet the plots manage to exclude as large 
a region in parameter space as in 0! As we discuss in Section V, the reason for this is the criterion used in |^ to 
identify a detectable signal. 

As seen in the bottom panel of Table ITTll the 7-ray emission towards the Galactic Center falls by at least two orders 
of magnitude on moving from a cusped NFW halo to a cored isothermal model. The point that the signal from the 
Galactic Center depends sensitively on the assumed halo profile, and so may have been overestimated, has also been 
made recently by Stoehr et al. j^. For example, an optimistic result for the 7-ray flux towards the Galactic Center 
was obtained in by using a cusped NFW model normalised to satisfy two constraints on the halo mass AI and 
circular speed Wh, namely 

M(r < 100 kpc) = (6.3 ±2.5) X 10"Mo, t>h(i?o) « 128 - 207kms-^ (10) 

If the normalisation is set to obtain the maximum flux (as is done for example in Fig. 9 of 0), then the models possess 
a local dark matter density substantially in excess of the usual value of ~ 0.3 GeV cm~^. Anyhow, even accepting 
the debatable proposition that the Galaxy did once have a pristine dark halo of the NFW form, the formation of 
the Galactic disk, bulge and bar will have substantially re-processed the dark matter distribution. Certainly, the 
evidence from the microlensing optical depth towards the Galactic Centre and the pattern speed of the Galactic bar 
are inconsistent with models dominated by dark matter in the central regions [T^ . 

Finally, let us illustrate how to convert the numbers in Table IIIII into photon fluxes. This requires adopting 
characteristic values for the particle physics parameters; here, we take — 100 GeV, and N^{av) = 10^^^ cm^^ s^^ 
for energies in excess of 1 GeV. Assuming the field of view is fixed at 3° (in other words, the semi- angle of the cone 
whose axis points towards the center of the objects is 1.5°), the implied photon fiuxes from the Sagittarius, Draco 
and Canis Major dSphs are as given in Table HVl 
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Typical Expected Flux Values (cm ^ s ^) for i5 > 1 GeV 



Model 


Galactic Center 


Sagittarius dSph 


Draco dSph 


Canis Major dSph 


NFW 


3.2 X 10"^ 


1.1 X 10""' 


9.9 X 10"^^ 


7.8 X 10"" 


Power-law 


6.0 X 10"" 


5.3 X 10"" 


5.1 X 10"^^ 


2.9 X 10"" 



TABLE IV: We show the annihilation 7-ray fluxes (cm s ^) for continuum emission above 1 GeV from the Galactic Centre 
and the three dSphs, assuming a field of view of 3°. These numbers are computed from the values of (J) given in Tableimiunder 
further assumptions of characteristic values for the neutralino mass and cross-section (viz. = 100 GeV and N^{av) = 10~^^ 
cm'^ s"^). In each case, results for a cusped NFW (7 — 1) and a cored power-law (a = 0) model are given. 



IV. COMPUTATION OF THE BACKGROUND 



There are three sources of background for the signal under consideration: hadronic, cosmic-ray electrons and 
diffuse 7-rays from astrophysical processes. The last is negligible for ACTs, but is the only one present for satellite 
experiments like GLAST or EGRET. Let us consider each source of background in turn. 



Hadronic and Electronic 



Bergstrom et. al. |3| use data taken with the Whipple ACT to derive the following expression for the hadronic 
background: 

'^^''^'^ {E > E,) = 6.1 X 10-3,,^, ( (11) 



where Chad is intended to take into account improved hadronic rejection expected in future ACTs, but is at present 
set to unity. 

Showers initiated by cosmic-ray electrons are indistinguishable from gamma rays. This contribution to the back- 
ground for ACTs is, according to 2j (who cite [s^ for this purpose): 

d$,- ^ o f E, "^-^ 



3 X 10-^ ("tI^ V ■ cm-V^sr-i. (12) 



Ail V 1 GeV J 



B. Diffuse Emission 

The diffuse 7-ray background is usually taken to be dominated by the Galactic j^^l or extragalactic contribution, 
depending on whether the target location is the Galactic Center or at higher latitudes (& > 10°). For example, a fit 
to the EGRET data [s^ at 1 GeV (dominated by the Galactic contribution) is given in 

*| = ):0-(^)-"c„-V.-GeV-, (13, 

where iVo(^, h) is a factor in the range 1-100, with higher values for the central regions of the Galaxy. In 0, only the 
extragalactic contribution from EGRET, estimated in ,3^], is considered: 

,^ / 77 \ -2.10±0.03 

d1>diff 00 , « O.N ,, 10-9 / ^0 \ „_-2„-l„,.-l.,r„^,-l 



,^ , = (7.32 ± 0.34) X 10"^ cm"^s"'sr"^MeV 

<mdE ^ ' V451MeVy 

f -p \ -2.10±0.03 

w 1.4 X 10-M — ^ ) cm-^s-^sr-^GeV-^ (14) 

y 1 GeV / 

So, the spectral indices of the Galactic and extragalactic contributions are about -2.7 and -2.1 respectively. 

However, the separation between the Galactic and extragalactic background is not clear. For example in |39l |. the 
case is made for a very low extragalactic background. Studying the region around the Galactic poles (h « 90°), it 
seems that, even there, most of the contribution is of Galactic origin. In particular, the main contribution is not 
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isotropic but correlated with known Galactic tracers. The EGRET collaboration concede that any simple model for 
the diffuse background is unlikely to work for all points in the sky and at all energies |40|. 

To be conservative, we normalize the flux to the EGRET data above 1 GeV and choose a spectral index of -2.1 
which is the worst case: 

d$diff .JE 

The emission above 1 GeV in the region of our interest can be downloaded from the EGRET website ■ The exact 
values for the diffuse emission are 6.7 x 10~^cm~^s~^sr~^ at the location of the Draco dSph {I — 86.4°, b = 34.7°), 
3.18 X lO^'^cm-^g-isr-i at the location of the Sagittarius dSph {I = 5.6°, b = -14.1°), and 1.2 x lO^^cm-^g-isr-i 
at the Galactic Center. 

The diffuse emission is the only background for satellite experiments. Its large variation with Galactic coordinates 
can make a weak source in Draco relatively brighter than strong emission from the Galactic Center, overwhelming the 
numbers in Table ITTll For ACTs, however, the hadronic and electronic backgrounds are much larger and independent 
of Galactic coordinates, so the hierarchy from Table Hill is retained. So, this raises the possibility that the Sagittarius 
dSph might have a higher signal-to- noise ratio with ACTs, but the Draco dSph is more clearly seen from satellites. 

V. THE DETECTORS 

A. Minimum Detectable Flux 

For the dSphs, the minimum detectable flux is determined using the prescription that, for an exposure of t 
seconds made with an instrument of effective area Agg and angular acceptance AJl, the significance of the detection 
must exceed 5a: 



> 5. (16) 



+ *bg 

Here, <I>^ denotes the neutralino annihilation flux in cm~^s~^sr~^, while $bg is the background flux. Any detector sees 
some photons from both dark matter annihilation and background, so the error in the measurement is cx y/^-y + $bg 
and not cx -\/$^. As pointed out in _4JJ , the use of the latter formula overestimates the significance of any detection. 

When studying the signal from discrete 7-lines, d^hg/dfl is the background flux falling under the annihilation line. 
If the background has a differential spectrum (P^^g/dildE = NoE~^ , and if the energy resolution of the instrument 
is <je/E, then the background under a line at energy Eq (i.e., in the interval [E^ — ge, Eq + as] containing 68% of 
the signal) is given by 

^ = ^V+^x, (.,/£;, 5), (17) 

with 

V {oeIE^ S) = ^-—^ \ (18) 

For ACTs, the background is the sum of three different power laws; for satellites, only the diffuse background is 
needed. 

In the literature, a number of different algorithms are used to define a detection of dark matter annihilation. Some 
authors additionally require a minimum number of detected photons, though this number is set somewhat arbitrarily. 
In 0, the minimum number is 25 for ACTs and 10 for satellite experiments; in it is 100 for ACTs. In 
no minimum number of detected photons is required, which allows the possibility of a high significance detection 
with a tiny number of received photons. In P|, a completely different strategy altogether is used: constraints on 
supersymmetric parameter space are found by requiring that Draco's flux be less than the least significant detection 
(the Large Magellanic Cloud at Aa 0]) above 1 GeV, resulting in a minimum flux for detection of 10~^cm~^s^^. In 
this way, the noise enters Ivaearly int"he expression. This explains why Tyler excludes a larg e region in mSUGRA 
parameter space from the non-detection of dSphs, despite the fact that the values of the integral (j?)) towards the 
dSphs are quite low. 
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B. Performance of the Detectors 



The detector characteristics of the different experiments are summarised in Table For definiteness, we use 
An = fO-^sr « 0.1° for the angular average when considering ACTs (appropriate for energies ~ 100 GeV) or GLAST 
(10 GeV), and Afi = 10~'^sr « 1° for EGRET (10 GeV). Also important is the observation time, which is chosen as 
t w 1 yr for satellites. For the next generation ACTs, assuming four telescopes, we use an observation time t w 100 h 
and an exposure A^s = 4 x 10* cm^. This seems reasonable, as CANGAROO 1451 an d the last phase of HESS (4^ 
will have four telescopes, while VERITAS will have as many as seven j^]. MAGIC jig uses a single 17 m mirror and 
has roughly the same performance as next-generation ACTs, but with a reduced threshold of 30 GeV. 



VI. RESULTS 



The following plots show the parts of the supersymmetric parameter space that can be probed through the detection 
of a 7-ray signal from neutralino annihilations. We typically show the region to which GLAST and a generic second 
generation ACT will be sensitive. The plots found in '2^ show the 7-ray flux in cm~^s~"'^ against photon energy. They 
are not appropriate for depicting the exclusion limits from observations of different parts of the sky because the flux 
changes and hence so do the points representing theoretical models. We prefer to use the type of plot presented in 
with Nry{av) (which depends exclusively on the particle physics model) versus (although other quantities could 
be used as well). 

From eq. H16|l . we write the condition for detection in a more convenient way for the plots: 



, , Anml 25 + J625 + 100A17Acffi$bs , . 



Here, we see that increasing the angular acceptance Afi can increase the signal to noise ratio. In fact, if both signal and 
background are constant, the significance increases (and the minimum value of Ny{av) that can be probed decreases) 
as V Afl. However, the signal is not constant as Table IIIII shows, and the angular acceptance that maximizes the 
significance does not necessarily coincide with the minimum angular resolution of the detector. So the optimal strategy 
is to scan between the minimum angular resolution and the maximum field of view, choosing the field for which the 
signal to noise ratio is maximised. This depends on the position in the sky and on the type of profile. For instance, 
the distant Draco looks like a point source and the maximum signal is for the smallest angle possible. For Sagittarius, 
the optimal angle is 0.4° for cored profiles and the smallest possible for cusped profiles. In order not to put too many 
lines in the plot, we have avoided drawing all the halo types and show only the extreme cases. 



A. Discrete Lines 



The annihilation of two ncutralinos gives rise to two photons with energy « m^. The region probed by the 
different experiments is shown in Fig.^ Also shown are ~ 1500 points in the mSUGRA par ameter space that comply 
with all the accelerator limits (including b 57, {g — 2)^^ and other accelerator limits [43 that are incorporated in 
DarkSusy). All the points, bar five, have spin- independent cross-section with protons or neutrons below 10~^ pb, 
thus compatible with limits set by the Edelweiss nuclear recoil detector, but not the disputed signal claimed by the 
DAMA experiment Also, the upward-going, neutrino-originated, muon showers have a flux of < 10^ km~^ yr""'^ 
(which according to Kurylov and Kamionkowski is the limit set by super-Kamiokande) . 

As the figure shows, the discrete annihilation line is very unlikely to be observed, even with the next generation 
instruments. It is just about detectable for the most promising targets under the most optimistic assumptions - 
the Sagittarius or the Canis Major dSph galaxies assuming a Moore profile and using next generation ACTs. Other 
possible models (such as NEW or cored profiles) and targets (such as the Galactic Center) are much less propitious 
still. For GLAST, only one line is shown - namely that for the Canis Major dSph, but even this lies above all physical 
mSUGRA models and so provides no constraints. In particular, monochromatic lines from the Galactic Center are 
not visible to GLAST. The difference between this work and that of is that the latter authors took a very high 
dark matter concentration in the center (the profile is just NFW, but the constant in front is set to ensure maximal 
flux given two weak constraints on the mass and the rotation curve). This causes the 7-ray flux in monochromatic 
lines from the Galactic Center as computed by Q to be over two orders of magnitude greater than the values obtained 
in this paper. 
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XX -> 11 




FIG. 1: Exclusion limits for the discrete line xx ~^ 77- For ^-U the experiments, only the most favorable cases are shown. The 
green, red and blue points correspond to mSUGRA models with QcDmh^ in the range 0.005-0.2, as discussed in Section IIIA. 
The red points satisfy the more stringent WMAP constraints 0.09 < flcDmh^ < 0.13. The exclusion limits for xX ~^ -^7 ^-re 
very similar and not shown here. 



B. Continuum Emission 



The continuum emission comes from hadronization and subsequent pion decay. The programme DarkSusy |3,1 | uses 
results from the PYTHIA code 51] to compute the photon multiplicity for each neutralino annihilation. Experimental 
sensitivities are shown in Fig. [21 for continuum emission above 1 GeV and 50 GeV. 

The continuum emission above 1 GeV can yield some constraints. Although we have computed the curves for four 
targets (Draco, Sagittarius, Canis Major and the Galactic Centre) and for the full range of models in Section II, 
we give only the most promising results in the figures. The Draco, Sagittarius and Canis Major dSphs may yield 
interesting constraints - but only if their dark halo profiles are strongly cusped (the Moore and the NFW profiles both 
rule out some supersymmetric models). Unlike the case of the Milky Way, cusped profiles are still possible for the 
dSphs. Notice, however, that substituting cored power-law models for NFW or Moore profiles causes the exclusion 
limit to move well above the supersymmetric parameter space of interest. For > 1 GeV, only curves for GLAST 
are drawn, as ACTs are insensitive at such low energies. 

Also shown in the upper panel of Fig.[5]is a line corresponding to the Milky Way observed at medium latitudes with 
the wide field of view of GLAST, as first suggested by Stoehr et al. j^. (This line lies almost exactly on top of the 
line for the Sagittarius dSph in the upper panel). Here, the Galaxy has been modelled with an isothermal power-law 
model, as opposed to the cusped models preferred by Stoehr et al. We agree therefore with the suggestion of Stoehr 
et al. that this is a promising target, as irrespective of whether the Galaxy is cusped or cored, there are always useful 
constraints on the supersymmetric parameters. Unfortunately, this attractive option is only available to GLAST and 
not for ACTs. 

The lower panel of Fig. |21 shows the prospects for detection of continuum emission above 50 GeV. For ACTs, the 
Canis dSph is the best target, though a detectable signal will again only be measured if the density profile is strongly 
cusped. For GLAST, the Galaxy at medium latitudes again leads to some constraints, though not as strong as when 
continuum emission above 50 GeV is studied. 
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FIG. 2: Exclusion limits for continuum 7-ray emission above 1 GeV (top) and 50 GeV (bottom). Only the most favorable 
cases are shown. For > 1 GeV, only curves for GLAST are drawn, as ACTs are insensitive at such low energies. Above 50 
GeV, curves are shown for both GLAST and second generation ACTs. The green, red and blue points correspond to mSUGRA 
models with Qcoyih^ in the range 0.005-0.2, as discussed in Section IIIA. The red points satisfy the more stringent WMAP 
constraints 0.09 < fJcDM^^ < 0.13. 



VII. CONCLUSIONS 



If the dark matter present in the Universe is composed at least in part by the hghtest supersymmetric particle, 
then this could manifest itself via 7-ray emission from pair annihilations. It is clearly important to estimate the 
likely magnitude of the neutralino annihilation signal. It is also important to identify the likely locations and spectral 
regimes in which the signal should be sought. This paper has provided new estimates of the signal towards the 
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Galactic Center and the nearby dwarf spheroidals using a variety of models. 

There have been a number of recent calculations predicting that the neutralino annihilation flux from the inner 
Galaxy will be detectable with forthcoming satellites like GLAST and with second generation atmospheric Cerenkov 
telescopes (ACTs) "5, These calculations assume that the cusped Navarro- Frenk- White (NFW) models for the 
Milky Way halo hold good. This assumption is in contradiction with a substantial body of astrophysical evidence 
about the inner Galaxy [H Hi I1,'t| . In any case, even if the Milky Way halo was originally of NFW form, the 
formation of the disk and bulge will have reprocessed the primordial dark matter distribution i6\. In contradiction 
with earlier results, we do not find the prospects of detecting the annihilation flux from the Galactic Center to be 
particularly promising. In particular, the 7-ray line coming from the the 77 and final states is not detectable 
either with second generation ACTs or with the GLAST satellite. We caution that many of the recent estimates of 
high flux are sensitively dependent on the assumptions made regarding the innermost structure of the dark halo. Even 
the best numerical simulations have difficulty in resolving structures on scales less than 1 kpc, and so the inner profile 
is always found by extrapolation. 

The high mass-to-light ratios of the Local Group dwarf spheroidals (dSphs) makes them attractive targets. Cusped 
profiles like NFW are not presently ruled out for dSphs like Sagittarius or Draco. It may be that the visible dwarf 
galaxy lies entirely within the central parts of a cusped dark matter halo. If so, then the optimum targets are the 
Sagittarius and Canis Major dSphs. The detection of monochromatic lines is still extremely difficult, but the GLAST 
satellite may detect excess continuum 7 ray emission. This is of course a less distinctive feature than a sharp line. In 
particular, if the Sagittarius or Canis Major dSphs have a strongly cusped dark halo profile [p ~ r~^-^ or p ~ t-^^), 
then some regions of supersymmetric parameter space can be ruled out. Again, however, this conclusion only holds 
good if the dark halo profile is cusped. Using a cored isothermal-like model for the dark halo, even the Sagittarius 
and Canis Major dSphs may be invisible to GLAST and second generation ACTs. 

Unlike Bergstrom et al. 0, we do not find the Galactic Center to be a promising location. Partly, this is because 
we believe that the Milky Way does not have a strongly cusped profile based on the available astrophysical evidence 
dllllll. Partly, this is because Bergstrom et al. chose a generous overall normalisation anyhow — they used the 
NFW model corresponding to the maximum flux which satisfies two weak constraints on the mass and the rotation 
curve. Accordingly, the local dark matter density is as high as ~ 0.6 GeV cm^'^ in their model. When the circular 
velocity curve of such a halo is combined with that for the disk and bar, then it necessarily violates the constraint on 
the Galactic rotation curve in the inner parts. One important caveat of our results — however — is that the possible 
effects of a central black hole are not included in our calculations. Here, we merely note that the observability of any 
expected signal depends on the manner in which the black hole grows |53l l54| . 

Stoehr et al. 0] have also recently emphasised that the 7-ray emission from the Galactic Center may have been 
overestimated by the use of too strongly cusped profiles. They suggest that the galaxy at moderate latitudes (|5| > 10°) 
may also be a good target for detecting the continuum emission (they do not study the line emission). This is not 
really an option for ACTs with their small field of view. However, it is an attractive possibility for GLAST, as the 
continuum emission is detectable irrespective of uncertainties in halo structure. For ACTs, the best targets remain 
the Sagittarius and Canis Major dSphs. 

Very recently, the Large Magellanic Cloud (LMC) has been suggested as another likely target 29]. Judging from 
[5^ . the average mass to light ratio of the LMC within 8.9 kpc is only ^ 3 (as opposed to ^ 100 for the compact 
dSphs). This is an upper limit to the central mass to light ratio. In other words, much as in the Milky Way, 
dark matter dominates the outer parts of the LMC and is responsible for the asymptotic fiatness of the ratio curve. 
However, the central parts of the LMC are dominated by the luminous bar and disk. The assumption that the dark 
halo dominates the gravitational potential everywhere is therefore not valid. Hence, the procedure used in of 
fitting the rotation curve to a NFW dark halo is fiawed. The gravitational potential of the gas and stellar disk and 
bar simply cannot be ignored in the central regions. 
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High energy 7-ray detectors 





HESS (I) 


VERITAS 


MAGIC 


EGRET 


GLAST 


Energy 


40 GoV-10 TcV 


50 GcV-10 TeV 


30 GoV-10 TcV 


20 McV-30 GcV 


20 McV-300 GcV 


as/E 


« 10% 


« 15% 


«20% 


< 10% 


« 5% > 10 GeV 


Aett (cm^) 


4 X 10® (> 100 GeV) 


4 X 10* (> 100 GeV) 


4 X 10® (> 100 GeV) 


1.5 X 10^ 


10* 


3>mm (cm"^s"^) 


8 X 10"^^(> 100 GeV) 


9 X 10"^^(> 100 GeV) 


« 10"" (> 100 GeV) 


10"^(> 100 MeV) 


3 X 10-®(> 100 MeV) 


Ang. res. (single 7) 


< 0.1° at 100 GeV 


< 0.1° at 100 GeV 


?s0.2° 


< 5.8° at 100 MeV 


2° at 100 MeV 
0.1° at 10 GeV 


Field of view 


4.3° - 5° 


3.5° 


«5° 


0.5 sr 


2.4 sr 



TABLE V: Performance of the gammarray detectors. Numbers quoted correspond to 5cr sensitivity after 100 hours of observation for ACTs and 1 yr for GLAST. 



